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In this paper by studying p-band fermionic system with nearest neighbor attractive interaction 
we find translation symmetry protected Z2 topological superfluid (TSF) that is characterized by 
a special fermion parity pattern at high symmetry points in momentum space k = (0,0), (0,7r), 
(tt, 0), (tt, tt). Such Z2 TSF supports the robust Majorana edge modes and a new type of low energy 
excitation - (supersymmetric) Z2 link-excitation. In the end we address a possible realization of 
such interacting p-band fermions with Z2 TSF. 
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Superconductivity /superfluidity is a paired state in 
many-body charged/neutral fermion systems introduced 
by Bardeen, Cooper and SchriefFer (BCS). In BCS's the- 
ory, the glue for the Cooper pair is phonon - the collective 
mode of atoms in solids. To describe such a quantum or- 
dered state, a local order parameter Ak = (ckC_k) is 
introduced that differs in different superconducting (SC) 
states (for example, s-wave, p-wave, d-wave, etc.). Re- 
cently, people found that SC states with the same lo- 
cal order parameter may have different topological prop- 
erties, that leads to the concept of "topological super- 
conductivity (TSC)"— . According to the characteriza- 
tion of "ten- fold way" of random matrix^i, there are 
three types of TSCs in two dimensions : D-type chiral 
TSC without time reversal symmetry {px ± ipy-wave^, 
s-wave with strong spin-orbital coupling^, ...), C-type 
chiral TSC without time reversal symmetry [dx ± idy- 
wave^) and Dlll-type Z2 TSC with time reversal sym- 
metry {{px +ipy) X {px - ipj;)-wave SG^). 

Recently, rapid advances in trapping and manipulating 
ultrocold atoms in optical lattice have made it possible to 
simulate topological superfluid (TSF) - a kind of "SC" in 
a neutral fermion system. For the paired states in optical 
lattices, the glues cannot be the phonons. Instead, peo- 
ple have more degrees of freedom to tune the interaction: 
the tunable direct attractive interaction by Fcshbach res- 
onance (see reviews ir^i), the effective attractive inter- 
action of composite fermions in Bose-Fermion mixture 
from background Bosons^, the anisotropic dipole-dipole 
interaction^. Because the non-Abelian anyon of D-type 
chiral TSF may be designed as a decoherence-free qubit 
and plays an important role in the realization of fault- 
tolerant topological quantum computatiortiSiiii, people 
pay more attention on its realization. For instance, for 
spinless (s-band) fermions in an optical lattice with the 
nearest neighbor attractive interaction, the ground state 
may be chiral Px ± ipy-wave TS C^^'^'^ . 

In optical lattices, besides the tunable interactions, 
people may have an additional degree of freedom to ob- 
tain the TSC by considering p-band fermions^^. In this 
paper we studied (spinless) p-band fermions with nearest 
neighbor attractive interaction on a square optical lattice 
and found that the ground state may he a. Px y- Py-weive 
SF - the Px-hand fermions are paired into px-wa.ve SF 
state and the Py-hand fermions are paired into pj,-wave 



SF state. Such a px x py-wave paired state is translation 
symmetry protected Z2 TSF against arbitrary transla- 
tion invariant perturbations. In addition, we found that 
this Z2 TSF supports the robust Majorana edge modes 
and there exists a new type of low energy excitations 
- (supersymmetric) Z2 link-excitations. Let us explain 
why there exist Z2 TSF beyond "ten- fold way" : if one 
enforces translation invariant on the system, there may 
exist sixteen symmetry protected TSCs/TSFs on a two 
dimensional (2D) lattice that are characterized by a spe- 
cial fermion parity pattern at four high symmetry points 
(k = (0,0), (0,7r), (^,0), (^,7r))i^. 

Mean field phase diagram and Z2 TSF: In this paper 
our starting point is the following effective model of spin- 
less p-band fermions with on-site and nearest neighbor 
attractive interactions on square lattice 

i i i 

where px-hand fermions only hop along x-direction and 
Pj,-band fermions only hop along y-dircction. U and 
Uq arc the on-site interaction strength and the near- 
est neighbor interaction strength, respectively [U > 0, 
Uq > 0). We denote orbitals index by pseudo-spin 
a = (px, Py) = (ti i) a-nd represents two sites on a 
nearest-neighbor link. In this paper we only consider the 
half filling case with ^ = 0. 

The model in Eq.([T]) is unstable against superfluid (SF) 
orders that arc described by and Ap for s-wave and 
p-wavc order parameters, where Ag = {ci-\-Cii), Ap = 
(cj-fCi+e^-t-) = (ciiCi+e^i), respectively. Then the effective 
Hamiltonian becomes 

■i i i 

+ t E ^y^+-yi - ^oAp E ^l4+e,i + h.C. (2) 
i i 

One can see that px-hand fermions only pair along x- 
direction, (ci-|-Ci+e^t) 7^ 0, {ci^Ci+ey\) = and py-h&nd 
fermions only pair along y-direction, (c^Ci+ej^j.) ^ 0, 
{ci^Ci+e.^l) = 0. For this reason we call the p-wave SF 
as Px X Py-wave state. In the SF state, the total fermion 
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FIG. 1: Phase diagram: region I and region II are trivial SC 
states with As 7^ 0, Ap = and As 7^ 0, Ap 7^ 0, respectively; 
Yellow regions (region III and region IV) are Z2 TSFs with 
As 7^ 0, Ap / and As = 0, Ap / 0, respectively; region V 
is a gapless state that we are not interest in. The blue lines 
with dots are the super-Z2 TSFs. 



number Np ~ "-i conserved only mod 2. Thus the 
fcrmion parity Pp = (—1)''^^ is good quantum number. 

After solving self-consistent equations by minimizing 
the ground state energy, we arrive at the mean-field phase 
diagram shown in Fig.l. In Fig.l, there are five regions 
: I, II, III, IV, V. To classify the topological proper- 
ties of the ground states in the five regions, we calculate 
fermionic parities at high symmetry points by defining 
four Z2 topological invariablesi^: for k=(0,0), {tt,tt), 
we have a trivial result as Ck=(o,o) = Ck=(7r.7r) = 0; for 
k=(0,7r), (tt, 0), we have 

Ck=(o,.) = Ck=(.,o) = 1 - e + - 4^2) (3) 

where 9(a:) = 1 if a: > and 9(x) = if x < 0. Thus we 
identify two distinct SF states : SF with trivial topolog- 
ical indices Ck=(o,7r) = Ck=(7r,o) = in region I (A^ 7^ 0, 
Ap = 0) and region II (A^ ^ 0, Ap 7^ 0) and Z2 TSFs 

with nontrivial topological indices Ck=(o, it) = Ck=(7r,o) — 1 
in region III (A^ 7^ 0, Ap 7^ 0) and IV (A^ = 0, Ap 7^ 0). 
From topological indices of Z2 TSF (OllO)r^ we found 
a special fermion parity pattern: even fermion parity at 
k = (0, 0) and k = (tt, tt) and odd fermion parity at 
k = (0, tt) and k = (tt, 0). This is the reason why we call 
it Z2 TSF. The SF states in region V (A^ ^ 0, Ap 7^ 0) 
is gapless which we are not interest in. In the following 
parts, we will focus on Z2 TSFs in region III and region 
IV. 

Majorana edge state: In this part, we study the edge 
states of Z2 TSFs. Fig. 2 illustrates the edge states in 
region III and IV. In the region IV (A^ = 0, Ap 7^ 0), 
Z2 TSF behaviors likes two coupled ID p-wave SCs, one 
along x-direction, and the other along y-direction. Be- 
cause each ID p-wave SCs with effective Hamiltonian 

HeS = t J2i ^\/U^t/U+e:./y - ^pC^OCt/;,»Ct/i,i+e,/^ + h.C. 

is really a Majorana chain, we get zero energy states at 




FIG. 2: Left figure: Edge states of TSF in region III which 
have a symmetry protected zero energy mode at k^: — 0, tt; 
Right figure: Edge states of TSF in region IV which have a 
fiat band. 

the boundaryi^. Indeed, the zero energy states can be 
mapped to a Majorana fermion system with a flat band 
exactly (See right figure in Fig. 2). In the region of III 
{As 7^ 0, Ap 7^ 0), the Majorana fermion system is still 
stable but has dispersionii (See left figure in Fig. 2). 

On the contrary, in the region I and II, there are no 
edge states at all. 

Z2 link- excitation: Next, we study the low energy ex- 
citations in the Z2 TSF states in region III and IV. One 
excitation is the Bogoliubov quasi-particle which has an 
energy gap Eqp in Z2 TSFs (See blue line in Fig. 3). An- 
other excitation is the collective mode - the Goldstone 
mode due to spontaneous global U(l) symmetry break- 
ing. In this part we emphasis that there may exist an- 
other type of low energy excitations for TSF on a square 
lattice " the Z2 link-excitation. 

A Z2 link-excitation is defined as the sign-switching 
of SC order parameter on a link position I = (/, J) = 

{io, io + ex/ey) as Ap { f^^' , .^'^ . From this def- 

Ap, 2 7= Zq 

inition, one can see that its particle number is con- 
served mod 2. This is the reason why we call it Z2 
link-excitation. By exact diagonalization calculations, 
we obtain the excited energy Eunk for a single Z2 link- 
excitation in Fig. 3 (See magenta line). From it one can 
see that the Z2 link-excitation even have lower energy 
than quasi-particle excitation and becomes relevant to 
low energy physics. 

We point out that an important property of Z2 link- 
excitation is the existence of a bound state (zero mode) 
on it for the Z2 TSFs states both in region III and IV. 
The inset of Fig. 3 shows the bound state on a Z2 link- 
excitation in the region IV. Then we define that is 
the energy for the empty bound state (or the even parity 
state) and E_ is the energy for the occupied bound state 
(or odd parity state). In the region IV, when UoAp 7^ t, 
there may exist energy splitting between E^ and due 
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FIG. 3: Three energy scales in the TSF states ; Eqp, Eunk, 
AE (The scale is t). The left inset is the scheme of the energy 
spectrum of quasi-particles, of which Eqp is the energy gap of 
Bogoliubov quasi-particle and AE is the energy splitting of 
the bound state of Z2 link-excitation. The right inset figure 
shows the bound state on a Z2 link-excitation with U/t = 
0.3, Uo/t — 3.25. Eiink is the excited energy of a Z2 link- 
excitation. 



to quantum tunneling effect, AE = E+ — By exact 
diagonalization calculations, wc derive AE in Fig. 3 (See 
dash violet line). The occupied bound state with odd 
fermion parity has lower energy than the empty bound 
state with even fermion parity, AE > (we have set 
the chemical potential to be zero). Now the Z2 link- 
excitation can be regarded as a fermion when the energy 
scale is smaller than AE. 

Particularly, in the region IV, when UoAp = t. the 
energy splitting of the two bound states becomes zero 
AE = 0. For this case, the energy dispersion of quasi- 
particles becomes a constant or the Bogoliubov quasi- 
particles have a flat band, that corresponds to a quasi- 
particle with divergent mass. Thus the quantum tunnel- 
ing effect by exchanging quasi-particles is frozen and the 
energy is degenerate for an occupied zero mode and an 
empty zero mode. Now the I2 link-excitation of even par- 
ity bound state is a Boson and that of odd parity state is 
a Fermion, that is really emergent super symmetry. And 
we call this TSF with supersymmetric I2 link-excitation 
to be super-'L2 TSF. 

We also study the Z2 link-excitations in region III and 
find that their properties are the same as those in re- 
gion IV (sec the dotted line in region III of Fig.l that 
corresponds to the super-Z2 TSF). 

Z2 vortex: After recognizing the properties of Z2 link- 
excitation, we turn to study Z2 vortex. For a SF on 
square lattice, Z2 vortex is really a 7r-flux on a plaque- 
tte which is confined at zero temperature and cannot be 
real excitation. On one hand, the Z2 link-excitations can 
be regarded as bound states of a pair of Z2 vortices on 
nearest neighbor plaquettes (A and B) (See Fig. 4). On 
the other hand, the Z2 vortex can be regarded as the end 



FIG. 4: Left figure: Illustration of the relationship between 
fermionic Z2 link- excitation (the green link) and e-type Z2 
vortex, m-type Z2 vortex (the shadow plaquettes) of Z2 TSF. 
There are two sub-plaquettes, one for e-type Z2 vortex, the 
other for m-type Z2 vortex; Right figure: Illustration of the 
relationship between supersymmetric Z2 link-excitation (the 
green link) and two-component super-Z2 vortex (the shadow 
plaquettes) of super-Z2 TSF. 



of non-local string of Z2 link-excitations and is generated 
by a string of link-excitations. 

When Z2 link-excitation is a fermion (UoAp ^ t), due 
to fermion parity constraint, Z2 vortex cannot hop from 
a plaquette to nearest neighbor plaquette by creating an 
additional fermionic Z2 link-excitation which violates the 
fermion parity conservation. Thus there must exist two 
types of Z2 vortices: e-type Z2 vortex on A sub-plaquette 
and the m-type Z2 vortex on B sub-plaquette. Each type 
of Z2 vortex lies on each sub-plaquette. In addition, there 
is mutual-semion statistics between e-type Z2 vortex and 
TO-type Z2 vortex. And an e-type Z2 vortex and an m- 
type Z2 vortex annihilate with each other into a fermionic 
Z2 link-excitation. Now the system has 4 superselection 
sectors: 1 (the vacuum), e (e-typc Z2 vortex), m (m- 
typc Z2 vortex) , tJj (fermionic Z2 link-excitation) and the 
fusion rule is e x e = m x m = ■0 ^ V' — 1, e x m = 
-0, e X ~ m, mx ijj — e. The fusion rule is the same 
as that of the Z2 topological order in toric code model 
and Wen's plaquette mode l^^i^^ . 

For the super-Z2 TSF, we find that Z2 vortex has inter- 
nal degree of freedom - it can be either e-type Z2 vortex or 
TO-type Z2 vortex. That is each Z2 vortex has two compo- 

nents as cr = ( ^ ). When same types of Z2-vortices (ee, 

toto) annihilate each other, we have a bosonic Z2 link- 
excitation which belongs to the vacuum sector; When 
different types of Z2-vortices (em) annihilate each other, 
we have a fermionic Z2 link-excitation. For this reason 
we call it super-Z2-vortex. See the schemes in Fig. 4. The 
nonAbelian representation of the fusion rule for the su- 
persymmetric case is a X a = 1 + ip where cr is a ma- 
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trix. This equation is much similar to that of the non- 
Abehan TSG^i^i where the vortex section a is not a 
matrix. Just for this reason, super-Z2 TSF may be an- 
other possible candidate for the fault-tolerant topological 
quantum computation. 

Physics realization: In general, to realize such ex- 
otic Z2 TSF, we need to design a spinless p-orbital 
fcrmionic model with significant nearest neighbor at- 
tractive interaction. Unfortunately, for p-band fermions 
the direct nearest neighbor attracting interaction is too 
small to be considered. A candidate is an interact- 
ing mixture of ultracold bosons on s-band and spin- 
less fermions on p-band in 2D square optical lattice^, 
of which the effective Hamiltonian is H ~ J2s=s pi^s — 

IJ'sNs) — Usp J2i f^s,inp.i where and Hp describes the 
hopping terms of s-obital bosons and p-obital fermions: 
Hs = -tsY.(i.j) ^I.^^sJ + ^ nsA^s^i - 1) and Hp = 

^W^i fi^fi+eA + ^11 J2i filJi+eyl ~ ^-L E j /jj., /i+Exi ~ 

*-L J2i fliJi+eyt' respectively, is the hopping pa- 
rameter for s-band bosons and is hopping pa- 
rameter for p-band fermions on neighboring sites par- 
allel/perpendicular to the bond direction. Pseudo-spin 
(t, i) denotes orbitals index {p^, Py). Due to t_\_ <^ t\\, 
we set t± to be zero. Usp > denotes the attractive 
interaction strength between s-band bosons and p-band 
fermions and Us > denotes the repulsive interaction 
strength between s-band bosons, fis is the chemical po- 
tential and Ns = ^ ng^i is the particle number operator 

i 

where fis^i = bl ^bsj and fip^i = fl^ji^ H- fl^Jti- Here, we 
consider the Mott insulator phase of s-band bosons with 
unit filling and a half filling case of p-band fermions. 

By tuning the interaction between s-orbital bosons and 
p-orbital fermions for the case of Us > Usp > 0, an 



s-orbital boson and a p-orbital fermion may bind into 
a composite p-orbital fermion with both effective on- 
site and effective nearest neighbor attractive interaction 
(See detailed calculations in Ref4). In this case, the 
corresponding parameters in Eq.[l] are t — > 2tst\^/Usp, 

u ^ Us - Wsp, Uo ^ -it^s/Us - ti/Usp - q/Usp, 

respectivelj^. Here is the operator of the composite 
fermion as bifi. From the fact Uq < 0, there always exists 
effective nearest neighbor attractive interaction that will 
leads to Px x py-wave pairing between composite fermion 
of the same p-orbital. If 4[/sp > Us > Usp, we have an 
additional on-site attractive interaction that will leads to 
s-wave pairing between composite fermion of the oppo- 
site p-orbital. When the effective nearest neighbor at- 
tracting interaction dominates, the ground state is just 
the Z2 TSF. 

Conclusion : In this paper we studied the interacting 
p-band fermions on square lattice and found the ground 
state may be translation symmetry protected Z2 topo- 
logical superfluid beyond the characterization of " ten- fold 
way" . Such Z2 TSF has even fermion parity at k = (0, 0) 
and k = (tTjTt), odd fermion parity at k = (0,7r) and 
k = (tt, 0). This Z2 TSF supports the robust Majorana 
edge modes and there exists a new type of low energy 
excitation - (supersymmetric) Z2 link-excitation. In the 
end we address a possible realization of such interact- 
ing p-band fermions - a mixture of ultracold bosons on 
s-band and spinless fermions on p-band in 2D square op- 
tical lattice. 
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